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ELASTIC  MODULI  OF  ALUMINUM-POLYXJRETHANE 
COMPOSITES  (ALUMERS) 


imoDucnoN 

In  acoustical  design,  it  is  often  desirable  to  have  an  acoustically 
transpairent  material  available;  i.e.,  a  material  with  density  p  (rho)  and 
dilatational  speed  c^  that  closely  match  the  density  and  compressional  wave 
speed  in  the  medium,  usually  (8ea)water.  This  ideal  may  be  closely  approached 
by  elastomers  in  the  rubbery  state.  These  are  sometimes  called  'rho>c  rubber” 
(a  trade  name  of  B.F.  Goodrich).  The  shear  modulus  G  of  these  elastomers  is  low 
and,  as  a  consequence,  they  do  not  offer  the  structural  support  that  may  be 
needed  in  some  applications. 

A  rigid,  acoustically  transparent  material,  based  on  a  fluoroepoxy  matrix 
into  which  glass  and  phenol  microballoons  are  dispersed,  has  been  described  in 
the  literature  [1] .  In  this  report  a  different  composite,  aimed  at  obtaining 
the  desired  match,  is  introduced;  it  consists  of  a  foamed  aluminum  matrix, 
impregnated  with  a  two-part  polyurethane  [2] .  The  term  "alumer”  (a  contraction 
of  aluminum  and  elastomer)  is  proposed  for  these  composites. 

It  may  be  noted  here  that  a  solid  material  necessarily  deviates  from  ideal 
acoustic  transparency  by  virtue  of  the  presence  of  a  finite  shear  modulus, 
causing  reflections  at  other  than  normal  incidence — even  with  perfectly  matched 
p  and  c^  [3]  .  This  effect  is  particularly  important  near  grasing  incidence  and 
near  thickness-shear  resonances. 

In  interpreting  the  experiments,  it  is  assumed  that,  at  low  frequency,  the 
composite  behaves  as  a  homogeneous,  isotropic  material,  characterised  by  two 
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independent  complex  moduli  (one  of  the  constituents,  the  foamed  aluminiim,  shows 
anisotropic  behavior;  it  is  assumed  that  this  does  not  make  a  noticeable 
difference  in  the  composite) .  One  may  compare  this  situation  with  the  case  of 
fluid-filled,  porous  materials  where  Biot*s  theory  predicts  the  presence  of 
three  waves:  two  in  the  solid  and  one  in  the  fluid.  The  coupling  between  the 
two  solids  encountered  in  this  study  demands  continuity  of  both  the  normal  and 
the  tangential  stress;  as  a  consequence,  it  would  be  difficult  to  imagine 
independent  waves  propagating  in  the  two  solid  constituents.  There  was  no 
literature  found  that  deals  with  this  possibility. 

This  report  discusses  the  measurement  of  the  bulk  modulus  of  pure 
polyurethanes  and  the  corresponding  alumers.  Possible  systematic  errors  in  the 
measurement  of  Young’s  modulus  are  discussed,  and  a  correction  model  is 
presented.  The  results  of  the  measurement  of  Young’s  modulus  E  of  the  alumers 
are  given  (the  Young’s  modulus  of  the  polyurethanes  is  taken  from  the 
literature) .  The  dilatational  wave  speed  is  computed,  and  the  resulting  pc- 
products  are  listed. 

BULK  MODULUS  OF  ELASTOMERS 

In  Fig.  1,  the  results  of  a  complex  bulk  modulus  K*  measurement  as  a 
function  of  frequency  are  shown  for  (soft)  natural  rubber  AA165-4  (see  Ref.  4 
for  a  discussion  of  the  meas\irement  method) .  The  modulus  is  expressed  in  terms 
of  real  and  imaginary  parts  by  K*  =  K'  *  iK' ' ;  the  loss  tangent  is  tan  = 
K' '  /K'  .  Only  the  solid  parts  of  the  curves  are  considered  to  represent  true 
values  for  K'  and  tan  One  sees  that  at  higher  frequencies  the  curve  for  K' 

(dashed  part)  turns  down;  this  is  due  to  the  fundamental  chamber  resonance  at 
about  8  kHz.  For  this  reason  and  because  of  another  structure  visible  in  most 
other  cases  (see  below) ,  the  validity  of  the  data  is  assumed  to  be  limited  to 
values  between  50  and  500  Hz.  The  pressure  in  the  measurement  coupler  is 
2.5  MPa  throughout. 
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FREQUENCY  (kHz) 


Fig.  1.  R«ault«  of  th«  dat»  raduetien  of  tk*  djnsmie  bnlk  Mdnln*  ■•aanreaant  ••  » 
fuaetioa  of  frequaaef  at  9*C  for  aatural  rabbor  AA.16S-4  (ogad  6  yoaro  la  ■aawatar)  .  Tba 
•olid  parta  of  tbo  earroa  (eoaaoetiag  tbo  aoaaarod  poiata)  aro  eoaaidarod  to  bo  tbo  real 
part  x'  aad  tbo  loot  taagoat  taa  6^  of  tbo  balk  aodalaat  tbo  daabod  parta  aro  affoetod  bp 
tbo  ebaabor  roaoaaaeo. 

One  expects  that  the  bulk  modulus  will  show  viscoelastic  behavior  in  its 
dependence  on  frequency  and  temperature,  although  less  so  than  the  Young’s  or 
shear  modulus.  In  Refs.  4  and  5  this  viscoelastic  behavior  is  analysed  for  some 
specific  elastomers  for  which  the  frequency  and  temperature  dependence  of  K' 
and  tan  showed  the  familiar  viscoelastic  structure.  In  the  present  study 
this  is  not  pursued  further;  instead,  an  average  is  given  in  the  tables  for  K' 
and  tan  over  the  frequency  range  of  50  to  500  Hz. 

As  an  indication  of  the  error  involved,  the  standard  deviation  of  the 
population  is  given,  estimated  from  the  variation  in  the  values  over  the  given 
frequency  range,  in  order  to  account  for  a  possible  systematic  frequency 
dependence;  a  large  value  of  this  error  is  generally  connected  with  a  trend 
visible  in  the  graphs.  For  mutual  comparison  of  the  experimental  K'  values, 
this  is  an  adequate  measure  of  uncertainty.  Since  the  method  relates  the  values 
of  K'  to  the  bulk  modulus  of  castor  oil,  the  absolute  value  of  K'  may  be 
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subject  to  a  systematic  error  of  about  2%  due  to  the  accuracy  limit  in  the 
compressibility  of  castor  oil  [8] . 

BULK  MODULUS  OF  PQLIUBEIHAMBS 

Figure  2  shows  the  values  for  K'  at  10,  20,  and  30* C  for  the  four 
polyurethanes  that  were  used  as  components  of  the  alumers;  FR  1526,  PR  1538, 

PR  1574,  and  PR  1500.  The  PR  designation  is  given  by  the  company  that  produces 
these  elastomers.  Polymer  Research  Co.  These  four  elastomers  were  chosen  based  on 
their  stiffness  at  room  temperature  and  at  a  frequency  near  1000  Hz;  PR  1526  is 
soft,  PR  1574  is  hard,  and  the  other  two  are  in  between.  A  description  of  their 
general  and  physical  properties  may  be  found  in  Ref.  7. 


Fis.  a.  BmI  pwt  of  the  bulk  aedulva  K*  averezwd  free  SS  to  SM  Ha  aa  a  fuaetiea  of 
taaparatura  for  four  polTvratbaaaa .  The  aolid  liaaa  eenaaet  tha  pointa  at  IS,  8S,  aad 
SS*0.  Tha  arror  bara  at  tba  aaaaaraaaat  poiata  rapraaaat  tba  atandard  dariatioa  of  tha 
population  (aatiaatad  froa  tba  rariatioa  with  fraquaaep)  at  that  ta^>aratura. 

The  solid  curves  were  drawn  to  guide  the  eye.  At  each  measurement  point, 
the  error  bar  represents  the  standard  deviation  for  the  population,  from  the 
average  over  the  frequency  range  froa  50  to  500  Hs.  One  can  see  a  similar 
dependence  on  temperature  for  each  of  the  four  curves.  The  values  of  the  bulk 
modulus  for  the  four  elastomers  at  a  given  temperature  vary  little,  no  more  than 
by  about  0.2  GPa,  as  is  typical  for  elastomers.  In  contrast,  the  values  for 
Young’s  modulus  vary  considerably  (see  the  relevant  section  below). 
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In  Fig. 3  the  values  for  tan  are  given  for  the  saae  polyurethanes  with 
the  same  paraseters.  The  error  limits  are  not  depicted,  since  the  overlap  would 
result  in  a  confusing  picture.  The  standard  deviation  in  tan  6^  is  typically 
0.004;  a  larger  value  generally  corresponds  to  a  noticeable  trend  in  the  data  as 
a  function  of  frequency. 
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?ig.  S.  Loaa  tangrat  tu  of  th«  bvlk  aodvlva  araragad  froa  6S  to  6SS  Ha  aa  a 
fonctioa  of  taapavatura  for  four  polTurathanaa .  thm  liaaa  eoaaaet  tba  poista  at  IS, 

2S,  and  3S^C. 

BUIX  MODULUS  OF  ALUMEBS 

Table  1  shows  the  results  for  the  values  of  K'  and  tan  at  three 

temperatures,  for  the  pure  polyurethanes  and  the  alumer  types  based  on  these 
elastomers  when  combined  with  foamed  aluminum  of  6X  nominal  density  and  pore 
sizes  of  10,  20,  and  40  pores  per  inch  (ppi) .  The  errors  given  are  estimates  of 
the  standard  deviations  for  the  population  obtained  by  averaging  over  the 
frequency  range. 
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Table  1.  Values  i'or  the  real  part  of  the  bulk  modulus  K'  and  the  loss  tangent 
tan  6^  for  pure  polyurethanes. 


Elastomer 

Pore  size  of 
foam,  ppi 

Temperature,  * 

C 

30 

20 

10 

PR  1526 

elastomer 

only 

K'. 

tan 

CP» 

2.23610.002 

0.00610.005 

2.33710.007 

0.01110.002 

2.48110.009 

0.01310.001 

10 

K'. 

tan 

GPa 

2.48610.007 

0.01110.002 

2.60110.010 

0.01610.002 

2.74710.010 

0.01810.002 

20 

K'. 

tan 

GPa 

2.52010.008 

0.01210.003 

2.65110.011 

0.01610.002 

2.84610.009 

0.01610.003 

40 

K'. 

tan 

GPa 

‘iL 

2.55610.006 

0.00710.003 

2.65410.008 

0.01110.002 

2.83610.007 

0.01410.003 

PR  1538 

elastomer 

only 

K', 

tan 

GPa 

2.28*10.005 

0.00410.004 

2.37110.003 

0.00710.003 

2.51110.004 

0.00910.003 

10 

K', 

tan 

GPa 

h 

2.54310.008 

0.00710.003 

2.66210.011 

O.OlllO.OOl 

2.71210.029 

0.03510.007 

20 

K'. 

tan 

GPa 

2.46810.010 

0.01010.004 

2.56810.017 

0.01710.003 

2.74410.017 

0.02210.003 

40 

K', 

tan 

GPa 

«K 

2.49110.007 

0.01110.002 

2.61010.005 

0.01210.003 

2.66210.031 

0.03610.002 

PR  1574 

elastomer 

only 

K', 

tan 

GPa 

2.44610.005 

0.00510.005 

2.53710.010 

0.00810.002 

2.68310.006 

0.00810.001 

10 

K'. 

tan 

GPa 

2.82410.011 

0.00610.002 

2.93810.016 

0.00410.002 

3.05310.020 

0.00510.004 

20 

K', 

tan 

GPa 

2.84210.018 

0.00310.003 

2.96210.017 

0.00310.001 

3.07610.019 

0.00310.002 

40 

K', 

tan 

GPa 

h 

2.85210.010 

0.00710.003 

2.97410.015 

0.01010.003 

3.11610.009 

0.01010.002 

PR  1590 

elastomer 

only 

K',  GPa 
tan 

2.36910.004 

0.00510.005 

2.46610.009 

0.01110.004 

2.61410.012 

0.01410.002 

10 

K',  GPa 
tan 

2.55510.00C 

0.01110.004 

2.71310.016 

0.01610.004 

2.84810.002 

0.02510.003 

20 

K',  GPa 
tan  d» 

2.60810.006 

0.00810.003 

2.75310.009 

0.01410.003 

2.87710.014 

0.02610.004 

40 

GPa 

2. 62110. X9 

2.73610.013 

2.91110.032 
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To  evaluate  the  influence  of  the  aluminuiD  matrix  on  the  effective  bulk 

modulus  in  the  alumers,  one  may  study  the  ratio  of  K  ' /Yi  where  the 

■  P  * 

numerator  is  the  bulk  modulus  of  the  pure  polymer  material  and  the  denominator 
the  bulk  modulus  of  the  given  alumer.  The  results  are  shown  in  Table  2.  In 
addition  to  the  results  at  a  given  temperature  and  pore  size,  averages  are 
shown  for  alumers  with  given  pore  size  and  different  temperatures  and  for 
different  porr  aizes  at  a  stated  temperature. 

The  results  in  Table  2  show  that  the  variations  with  pore  size  and 
temperature  of  K'  for  a  given  filler  appear  generally  insignificant,  except 
for  the  alumers  based  on  PR  1538  with  10  ppi  2uid  40  ppi  foamed  aluminum  at  10 
"C,  for  which  the  K'  is  significantly  larger  than  at  the  two  higher 
temperatures.  If  one  looks  at  the  average  over  all  the  data  for  an  alumer 
based  on  a  given  polymer,  one  sees  a  gradual  incre2ise  of  the  ratios  in  the 
following  order,  from  low  to  high:  PR  1574,  0.864;  PR  1526,  0.886;  PR  1500, 
0.008;  PR  1538,  0.016. 


One  may  give  the  following  interpretation  of  the  data  in  Table  2,  when 

averaged  for  each  of  the  four  groups.  Ignoring  the  (small)  compliance  of  the 

aluminum  in  the  matrix,  one  could  approximate  the  compliance  of  the  composite 

by  that  of  the  polymer,  assumed  not  to  be  influenced  by  the  matrix,  equal  to 

the  porosity  of  the  foamed  aluminum  f  multiplied  by  the  inverse  bulk  modulus 

of  the  elastomer.  Thus  K'  /K*  s  (V  P  )/\P  (V  +V  )]  =  4,  where  V  is  the 

p'  ^  '  P  *  P  '  ^  p'  ^  '  p 

volume  of  the  polymer  in  the  composite,  V  the  volume  of  the  foamed  aluminum, 

and  p^  the  specific  compliance  of  the  pure  polymer. 


As  a  consequence,  one  would  expect  that  the  above  ratio  of  moduli  would 
be  approximately  equal  to  the  porosity  of  the  foamed  aluminum,  nominally  equal 
to  0.94.  The  fact  that  the  ratios  are  smaller  than  this  number  indicates  that 
the  aluminum  matrix  decreases  the  compliance  of  the  polymer  filler. 
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Table  2.  Ratio  of  the  bulk  modulus  lor  the  pure  polyurethane  and  the  bulk 
modulus  of  the  altimer  based  on  that  polyurethane. 


Elastomer 

Pore  size  of 

Temperature, 

•c 

type 

foam,  ppi 

30 

20 

10 

Average, 

Average , 

pore  size 

alumer 

PR  1526 

10 

0.900 

0.899 

0.903 

0.901 

20 

0.887 

0.882 

0.872 

0.880 

40 

0.875 

0.881 

0.875 

0.877 

Average  for 
temperature 

0.887 

0.887 

0.883 

0.886 

PR  1538 

10 

0.897 

0.891 

0.926 

0.905 

20 

0.925 

0.923 

0.915 

0.921 

40 

0.916 

0.908 

0.943 

0.922 

Average  for 
temperature 

0.913 

0.907 

0.928 

0.916 

PR  1574 

10 

0.866 

0.863 

0.879 

0.869 

20 

0.860 

0.856 

0.872 

0.865 

40 

0.858 

0.853 

0.861 

0.857 

Average  for 
temperature 

0.861 

0.857 

0.871 

0.864 

PR  1590 

10 

0.927 

0.910 

0.918 

0.918 

20 

0.908 

0.896 

0.909 

0.904 

40 

0.904 

0.901 

0.898 

0.901 

Average  for 
temperature 

0.913 

0.902 

0.908 

0.908 
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In  the  literature,  one  finds  models  that  are  designed  to  relate  the  bulk 
and  shear  modulus  of  a  composite  to  the  properties  of  the  components.  The 
model  of  Van  der  Poel  [8]  (also  found  in  Christensen  [9])  mas  applied  here. 
The  values  found  are  considerably  farther  removed  from  the  measured  values 
than  those  that  follow  from  the  compliance  of  the  polymer,  as  given  above.  For 
instance,  in  the  case  of  the  alumer  consisting  of  6X  relative  density  aluminum 
and  PR  1574  filler  at  30*G,  the  value  of  K' ,  according  to  Van  der  Poel’s 
model,  is  3.9  GPa,  about  SOX  larger  than  the  experimental  value.  The 
experimental  values  for  K'  show  no  clear  correlation  with  the  shear  modulus; 
this  mould  appear  to  confirm  the  result  of  Van  der  Poel’s  model,  namely,  that 
the  effective  K'  depends  only  on  the  shear  modulus  of  the  matrix,  and  not  on 
that  of  the  filler.  None  of  the  models  foxind  in  the  literature  account  for  the 
pore  sise  of  the  matrix,  only  for  its  porosity.  Temperature  does  not  enter 
explicitly,  since  the  model  uses  the  values  of  K'  of  the  tmo  constituents  at  a 
given  temperature.  The  models  are  developed  for  the  steady  state,  but  the 
results  for  the  bulk  modulus  should  not  be  noticeably  different  from  those  at 
lorn  frequency. 

The  values  for  the  loss  tangent  tan  (Table  1)  shorn  some  systematic 
effects.  The  standard  deviation  of  the  population  varies  from  about  0.001  to 
0.004.  The  larger  values  usually  reflect  a  trend  that  is  visible  in  the 
graphs . 

For  the  pure  polyurethanes  at  10  and  20*  C,  the  values  of  tan  for  PR- 
1590  and  PR-1526  are  close,  so  are  the  other  tmo  (see  Fig.  3).  At  30*C  all 
four  moduli  come  close  together.  In  almost  all  cases,  tan  of  the  composite 
is  higher  than  that  of  the  elastomer  by  not  more  thaui  a  factor  of  about  2.  An 
exception  is  found  in  the  alumers  based  on  PR  1538;  at  10*G  the  loss  tangent 
is  up  to  a  factor  of  four  larger  for  the  alumers  than  for  the  pure  elastomer 
at  the  three  pore  sizes. 
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STRUCTURE  IN  THE  GRAPHS  FOE  BULK  MODULUS  AND  LOSS  TANGENT 

Figure  1  is  an  example  of  a  favorable  case  in  that  the  smooth  appearance 
of  the  curve  shows  the  variation  of  modulus  and  loss  tangent  with  frequency 
below  the  chamber  resonance.  For  the  majority  of  other  polymers  a  structure 
appears  in  the  K'  and  tan  curves,  as  illustrated  by  Fig.  4,  which  shows 
these  quantities  for  a  sample  of  PR  1538  at  three  temperatures.  The  structure 
is  even  more  pronounced  for  a  hard  elastomer  like  PR  1574,  shown  in  Fig.  5, 
and  it  is  found  for  all  the  alumers  that  have  been  investigated. 


0.4 

0.3 


0.2 

0.1 


0 

0.01  0.1  I.O  10. 

FREQUENCY  (kHz) 

Fig.  4.  Saaalta  of  tho  dot*  roduetioB  of  tho  balk  aodulao  aoaoaroMBt  for  PR  16SS  bo  b 
fuBctioB  of  froquoBej  Bt  tbroo  toapBrBtaroo .  Tb«  dBtB  Bbor*  SSS  Ho  btb  Bot  eoBoidorod  to 
roprooBBt  tbo  d/BBaie  balk  Bodulao,  oIbco  thop  or*  Bffoetod  bj  coBoaqaBBCBo  of  tbo 
ozporiaoBtBl  aothod.  Solid  cotto,  IS  *C|  doobod  cottb,  3S*0)  dottad-dBobad  eurra,  SS*C. 
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HKOUCNCy  (kHz) 

Fig.  8.  Saanlt*  of  tli*  data  radaetion  of  tha  bulk  aodalaa  MaanraMat  for  Pk  1874  aa  a 
fonetioB  of  fraqvaney  at  20*0.  Tba  data  aboTo  800  Ba  ara  not  eoaaidarad  to  rapraaant  tba 
djaaaie  balk  aodulaa,  alaea  tbap  ara  affaetad  bp  eoaaaqnaaeaa  of  tba  azpariaaatal 
aatbod. 


At  one  time,  it  was  thought  that  this  structure  sight  be  related  to 
relaxation  in  the  material  at  hand.  The  following  arguments  are  contrary  to 
this  notion: 

1.  The  width  of  the  loss  peak  is  much  narrower  than  would 
correspond  to  a  single  relaxation  time,  while  the  visco¬ 
elastic  phenomena  are  always  related  to  a  broad  spectrum  of 
relaxation  times. 

2.  The  structure  shifts  to  lower  frequencies  when  the  temperature 
is  decreased;  this  is  opposite  to  the  temperature- frequency 
shift  known  in  viscoelastic  behavior. 

The  most  probable  cause  for  this  structure  is  the  presence  of  a  bar-type 
resonance  in  the  (cylindrical)  sample.  Under  this  assumption,  the  resonance 
frequency  would  scale  according  to  Cg/C-  ~  V(E/p)/2.,  where  c^  is  the  bar  speed, 
t  the  length  of  the  cylinder,  and  p  the  density  of  the  material.  It  was  found 
for  a  sample  of  PR  1574  that  the  resonance  frequency  shifted  upward  by  a 
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factor  Tory  close  to  1.5  when  the  cylinder  length  sas  reduced  to  2/3  of  its 
original  length,  confiraing  the  assuaed  length  dependence.  The  structure 
disappeared  coapletely  out  of  the  giren  frequency  range  when  the  cylinder  was 
replaced  by  a  sphere  equal  in  diaaeter  to  the  cylinder.  On  the  other  hand,  the 
shift  of  the  resonance  frequency  did  not  appear  to  follow  the  VE  dependence. 
Despite  this  failure  to  follow  the  VE  dependence,  the  basic  cause  is  still 
thought  to  be  a  bar  resonance.  One  can  easily  advance  arguments  why  deviations 
froa  a  siaple  bar  resonance  would  occur:  the  cylinder  is  loaded  by  the  castor 
oil  fill-fluid  and  the  space  between  the  cylinder  and  the  chaaber  is  quite 
narrow,  such  that  viscosity  and  flow  about  the  cylinder  aight  play  a  role. 

In  suaaary,  these  arguments  lead  to  the  conclusion  that  the  observed 
structure  is  not  due  to  an  effect  of  molecular  processes  on  the  bulk  modulus. 
In  order  to  avoid  interference  of  the  observed  feature  with  the  stated 
aeasureaent  values,  the  averaging  of  K*  and  tan  5^  described  before  was 
limited  to  the  range  of  50  to  500  Hs,  where  this  type  of  structure  was  not 
visible  in  any  of  the  samples  investigated.  Some  remaining  deviations  from  a 
smooth  curve,  such  as  the  single-measurement  point  at  0.1  kHs  in  Fig.  4,  are 
considered  to  be  spurious. 

ELASTIC  FRQFEBIIES  OF  FOAMED  ALUMINUM 

The  measurement  of  Young’s  modulus  and  Poisson’s  ratio  of  various  types 
of  foamed  aluminum  is  described  in  Ref.  10.  The  relevant  properties  are  given 
here  (Table  3)  in  order  to  assess  the  way  in  which  these  properties  enter  into 
the  composite  alumer.  The  Young’s  modulus  measurement  consists  of  observing 
the  complex  ratio  of  the  accelerations,  or  velocities,  of  the  base  and  the  tip 
of  a  bar  as  a  function  of  frequency,  and  finding  the  dimensionless  wave  number 
for  extensional  waves  by  means  of  a  complex  root  finder  [11] .  In  the  study 
reported  in  Ref.  10,  the  velocities  were  measured  by  means  of  a  laser  Doppler 
vibrometer;  for  the  measurement  of  the  Young’s  modulus  of  the  alumers, 
discussed  below,  accelerometers  were  used  to  determine  the  accelerations.  The 
latter  arrangement  is  similar  to  the  one  described  in  Ref.  7,  which  was  used 
to  determine  Young’s  modulus  of  the  pure  polyurethanes. 
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Table  3  -  Elastic  properties  of  foamed  aluminum;  ppi  -  pores  per  inch; 

E  -  Young’s  modulus;  V  -  Poisson’s  ratio;  p  -  density;  ^  - 
porosity;  K  -  bulk  modulus. 


ppi 

E  (MPa) 

V 

P  3 
kg/m^ 

* 

K  (MPa) 

Dilat. 

modulus 

(MPa) 

Dilat . 

speed 

(m/s) 

Nominal  relatire  density  6X. 

10 

l/sjj  =  207 

0.33  = 

177 

0.03 

- 

— 

■“13/®!! 

20 

I/.33  =  1»4 

0.42  = 
"®13/®33 

160 

0.04 

120 

Cjj=104 

1053 

1/sjj  =  157 

0.26  = 

175 

0.035 

1100 

"*13/®11 

0.23  = 

-•l2/‘U 

40 

174 

0.332 

161 

0.04 

173 

260 

1271 

Nominal  relatire  density  12X. 

10 

1/Sjj  =  238 

0.24  s 

258 

0.004 

- 

- 

— 

“®13/®11 

20 

CO 

II 

B 

0.10  s 

275 

0.808 

- 

— 

- 

-®13/*11 

0.27  = 

-®12/®11 

40 

251 

0.22 

281 

0.806 

140 

287 

1011 
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In  Ref.  10,  it  was  mentioned  that  some  of  the  foamed  aluminum  samples  did 
not  appear  to  be  isotropic  to  visual  inspection.  An  attempt  was  made  to 
explain  the  measurements  in  terms  of  the  set  of  compliances  that  correspond  to 
uniaxial  symmetry.  Stiffness  and  compliance  matrices  are  given  in  Appendix  A. 
Where  applicable.  Young’s  moduli  and  Poisson’s  ratios  in  Table  3  are  expressed 
in  terns  of  the  compliances  s^^ . .  Only  in  some  of  the  cases  was  there 
sufficient  information  to  compute  the  dilatational  speed  in  the  foamed 
aluminum. 


The  values  for  Young’s  modulus  of  foamed  aluminum  reported  in  Ref.  10  may 
be  compared  with  a  quadratic  dependence  of  the  modulus  on  the  relative  density 
of  the  foam,  as  described  by  Friis  et  al.[12]  and  Warren  and  Kraynik  [13] 
(these  authors  also  give  theoretical  expressions  for  Poisson’s  ratio).  Warren 
and  Kraynik  present  a  justification  for  this  functional  dependence  by  a  model 
based  on  a  tetrahedral  unit  cell,  which  contains  four  identical  half-struts 
joined  at  equal  angles.  For  low-density  foams  the  Young’s  modulus  £,  relative 
to  the  modulus  of  the  solid  material  E^,  is  given  by 


(1) 


where  is  the  relative  density  and  C  is  a  constant  that  depends  on  the 
assumed  strut  geometry:  C  =  0.91  for  circular  cross  section,  1.10  for 
triangular  cross  section,  and  1.53  for  a  Plateau  border  (the  space  between 
three  identical,  mutually  tangent  circles,  typical  for  liquid  foams  of  very 
thin  films) . 


Figure  6  shows  the  values  for  E  given  in  Table  3  as  a  function  of 
relative  density  p^  with  three  curves  according  to  Eq. (1) ,  with  C  =  0.667, 
1.00,  and  1.35.  The  value  of  0.67  for  C  gives  reasonable  agreement  between  the 
theoretical  curve  and  the  measurement  points,  at  least  for  the  points  of 
nominal  relative  density  of  6X. 
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Tig.  6.  UMawad  taIm*  of  Tatwg'a  aodslM  B  for  fo*a*d  alwalavst  ralatlra  to  tk* 
■odulua  pf  th«  aolid  aatpripl  B^,  pp  p  foBetion  of  rplptirp  dpBsity  p^,  eeapprpd  witk  p 
qppdrptie  Bodpl  B/B  m  0  p  ^  (Wprrap  pad  IrPTaik  [XI])  |  i^ptp  p  la  tka  ralptirp  daaaitj 
pad  0  p  ceaatpat  thpt  dapaada  oa  tha  aaaaaad  atraad  eroaa  aaetioa.  Solid  liaa,  C  >  X.| 
daaltad  liaai  0  ■  l.SS|  dettad-daa^ad  liaa,  0  *  S.ST.  Beaiapl  ralatira  daaaitpi 
triaaglaa,  eS|  aquaraai  ISB. 

Warren  and  Kraynik  discuss  another  aodel  based  on  a  cubic  unit  cell.  The 
resulting  dependence  of  E  is  linear  in  p^.  This  fits  the  slope  of  the 
experinental  points  better,  but  the  cuaerical  ralues  are  off  by  an  order  of 
aagnitude . 

T01]N6*S  MODULUS  FOE  FUSE  PQLTUBEIBANES 

Figure  7  shows  the  real  part  E'  of  Young’s  aodulus  as  a  function  of 
frequency  for  the  four  polyurethanes  used  in  this  study  at  22*0  teaperature. 
The  loss  tangent  of  these  polyaers  is  giren  in  Fig.  8.  These  graphs  were 
coaputed  froa  the  polynoaial  coefficients  giren  by  Capps  [7] .  One  sees  that 
PR  1526  is  soft  and  has  the  largest  loss  tangent,  PR  1538  and  PR  1500  hare 
coaparable  E' ,  and  PR  1574  is  considerably  harder. 
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Fig.  7.  E«»l  p»rt  of  Touas'a  Mdulua  S'  m  a  fnaetioa  of  fraqMBCjr  at  93*C  for  four 


polTwathaaaa ,  ooi^tad  from  polxiinaial  eoaf f lelaata  ia  Oappa  [7]  . 
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SYSIEMAXIC  gPMM  IN  TOQNG’S  IttOXHAJS  MBASOIBffiNT 

The  measurenent  of  Young’s  modulus  of  the  mlumers  was  performed  without 
the  benefit  of  an  enrironmental  chamber,  at  room  temperature  which,  of  course, 
fluctuates  and  is  a  source  of  systematic  error  in  the  data. 

By  deliberately  changing  the  temperature  of  the  room,  it  was  possible  to 
obtain  an  estiaiate  of  the  effect  of  the  temperature  fluctuations.  The  ralue  of 
E'  decreases  for  increasing  temperature,  concurrent  with  riscoelastic 
behavior,  but  the  shift  is  larger  than  expected  on  the  basis  of  the 
temperature  effect  computed  from  the  information  in  Ref.  7.  The  change  in  the 
loss  tangent  is  comparable  to  that  of  the  pure  elastomer  according  to  Ref.  7. 

In  the  listing  of  the  data  below  it  will  be  assumed  that  the  nominal 
temperature  is  22*G.  The  temperature  fluctuations  are  estimated  to  result  in 
an  uncertainty  in  £'  of  t6  UPa  and  ±0.006  in  tan  ffg. 

In  almost  all  of  the  results  for  the  real  part  and  loss  tangent  of 
Young’s  modulus,  one  notices  a  structure  at  regular  intervals.  A  typical 
example  is  found  in  the  reduced  data  for  the  pure  polyurethane  PR  1574,  Figs. 
9  and  10.  While  investigating  this  effect,  it  was  found  that  other  modes  of 
bar  vibration,  like  flexural  and  torsional  modes,  did  not  appear  to  cause 
these  measurement  problems.  The  features  shown  in  Figs.  9  and  10  appear  near 
values  of  k£.  (k  is  the  complex  wavenumber) ,  for  which  the  absolute  ratio  of 
the  accelerations  of  base  and  tip  is  maximum.  A  possible  cause  for  this  effect 
is  derived  as  follows. 
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Fig.  fl.  Baaolta  of  tk«  dst*  r«d«etiaa  for  tk«  m*1  part  of  Tevag'a  aednlaa  B*  of  FI  1B74 
polparatlMUM  at  $3*0  as  a  faaetloa  of  fratiomejr.  TIm  dUrriatloas  fro*  a  saootht  BOBetoaie 
earr*  mrm  mmmvmmd  to  bo  doo  to  oxporiMatal  offoets  dlseossad  ia  tbo  toxt. 


FREOUENCV  (kHi) 

Fig.  Ig,  Booolto  of  tho  data  radoetioa  for  tbo  loos  taagoat  taa  of  Touag's  aodalao  of 
FB  1C74  poljrarotbaao  as  a  faaetioa  of  frogaoaep.  Tbo  doriatioas  from  a  soooth  earro  aro 
aosaoad  to  bo  dao  to  oaporiaoatal  offoeto  dioeaoood  ia  tbo  toxt. 
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The  observation  equation  for  the  reduction  of  the  data  is  given  bj 

R  =  cos  s  -  r  B  sin  b,  (2) 

where  R  is  the  complex  ratio  of  the  acceleration  (velocity  or  displacement)  of 
the  base  to  that  of  the  tip  of  the  bar;  and  b  is  the  complex  dimensionless 
wavenumber,  b  =  \JL.  The  symbol  r  represents  the  ratio  of  the  masses  of  tip  and 
bar,  T  =  approximating  the  errors  in  b  and  R,  Ab  and  AR,  by  the 

first  term  in  a  Taylor  series  expansion,  one  can  relate  the  relative  errors  by 

Ab  ^  _ cos  B  -  rB  sin  b _  ^ 

B  "  ■  .  V  2  R. 

B  Sin  B  (1+r)  •«-  r  B  cos  b 

This  expression  shows  that  near  a  maximum  of  R  the  numerator  is  finite, 
while  the  denominator  goes  to  sero.  (To  simplify  the  argument  it  is  assumed 
that  the  loss  tangent  is  sero.  For  complex  values  of  R  and  b  with  small  loss 
tangents,  the  thrust  of  the  derivation  remains  valid)  .  Thus  near  maxima  of  R 
one  finds  a  considerable  amplification  factor  resulting  in  a  large  relative 
error  for  k£.  given  a  moderate  relative  error  in  R  (see  Frits  [14]  ;  Capps  [7] 
gives  references  to  other  error  analyses).  On  the  other  hand,  for  values  of  s 
where  R  is  close  to  sero,  the  resulting  error  in  b  goes  to  sero.  For  a  small 
loss  factor  the  relative  error  in  E'  is  the  negative  of  twice  the  relative 
error  in  s,  and  thus  for  complex  R  and  s,  one  has 


AE  cos  B  -  rs  sin  s _  „  AR 

2  ^ 

E  (l-»-r)s  sin  s  -i-  r  s  cos  b  R. 

Analysing  the  expression  on  the  left  side  of  this  equation  in  terms  of 
the  real  part  of  Young’s  modulus  E*  and  loss  tangent  tan  one  finds  that 
for  a  small  loss  factor  the  relative  error  in  E'  and  the  error  in  tan  A«  are 

a 

approximately  given  by 


and 


I 


A (tan  Ag) 


(5) 

(6) 
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where  Re  and  la  indicate  the  real  and  imaginary  parts  of  a  complex  quantity. 

One  notices  also  that  for  a  0  the  error  in  E'  increases  beyond  bounds. 
This  causes  the  random  error  to  fluctuate  more  extensirely  at  low  frequency; 
in  some  cases  where  the  random  error  is  small,  one  finds  a  corresponding 
systematic  deriation  in  E' . 

Figures  11  and  12  show  the  computed  beharior  of  the  errors  in  E'  and  tan 
dg  ^or  a  IX  error  in  R,  as  a  function  of  frequency,  for  the  above  example  of 
pure  PR  1574  polyurethane.  The  errors  have  the  location  and  appearance  of  the 
structure  in  the  reduced  data  from  the  experiment  (Figs.  0  and  10). 


nCQUENOr  (MU) 


PiS>  Oeapvted  pweaBtag*  error  As' /■'  la  tbe  real  part  of  Teaag'o  aedalao  ■'  for  a 
18  error  ia  the  ratio  8  of  the  aeeeleratioao  of  baoe  aad  tip  of  the  bar  for  polpurethaao 
1S74 
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nKoucNor  (utt) 


Pig.  i3.  Coap«t«d  •rrer  A  (tu  A  )  is  ths  leas  tssgsat  tss  5  of  Tesag's  sedslss  for  a 
iX  error  is  the  retie  I  ef  the  eeeeleretieas  ef  base  aad  tip  of  the  bar  for  polporathaso 
FI  1574 

Figure  13  shoirs  the  effect  on  tan  6^  of  correcting  the  data  near  the 
first  occurring  structure  by  applying  percentage  changes  in  R.  A  value  of 
-1.4X  appears  to  reduce  the  systematic  deviation  to  a  random  fluctuation.  Even 
a  small  variation  of  the  relative  error  produces  a  noticeable  deviation.  The 
corresponding  effect  of  this  -1.4X  correction  in  R  on  S'  is  shown  in  Fig.  14. 
This  same  correction  of  -1.4X  carried  through  the  higher  frequency  range  does 
not  give  smoother  curves  for  £'  nor  for  tan  for  the  structures  at  higher 
frequency. 
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FNCOUENOr  (kHi) 

fig.  ii.  of  opplfiac  Toriowo  poreoatos*  ekaagao  to  tha  ratio  1  oa  tha  atractara 

ia  taa  aaar  1.6  kla  ia  Pig.  18 


FRCOUCNCr  (hHz; 


Pig.  14.  If  fact  of  applfiag  rarioaa  pareaataga  ekaagao  to  tha  ratio  ■  oa  tka  atraetara 
ia  ■'  aaar  1.6  Ua  ia  Pig.  6 
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In  some  cases,  a  systematic  Tariation  of  E'  at  the  low  frequency  end  of 
the  range  was  apparent,  amenable  to  correction  by  means  of  a  percentage  change 
in  the  ratio  R.  Figure  15  shows  an  example  for  the  alumer  consisting  of 
PR  1528  in  a  matrix  of  6X,  40  ppi  foamed  aluminum.  It  appears  that  AR/R  =  -2% 
girea  the  best  result.  This  same  percentage  correction  was  applied  throughout 
the  whole  frequency  range,  with  some  success  at  the  first  peak,  but  with 
little  effect  on  the  ones  at  higher  frequency.  The  effect  on  tan  is  quite 
modest:  the  low  frequency  range  is  hardly  affected,  the  first  peak  in  the  loss 
tangent  data  between  1  and  2  kHs  appears  to  have  been  reduced,  but  the 
following  peaks  are  again  practically  unaffected  by  this  correction. 


Fig.  IS.  of  oppl/ing  m  eorroetlon  of  -2%  to  tho  ratio  K  oa  tho  ralaao  of  B'  (from 
tho  data  raduetion)  ao  a  faaetioa  of  fraqaaaep.  Tha  oaapla  io  PB  1538  in  eX,  48  ppi 
foaaad  alaaiaan.  Solid  earra,  aaeorraetad|  daokad  earra,  eorraetad. 

Before  discussing  possible  causes  for  an  error  in  R,  the  consequence  of 
an  error  in  the  ratio  r  =  explored.  Again,  approximating  the  errors 
by  the  first  term  in  a  series  expansion  one  finds  that 
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^  _ _ r  sin  « _  ^ 

B  (1+r)  sin  a  ■»•  r  s  cos  a  r,  ' 

where  Ar  is  the  error  in  r. 

Bj  a  procedure  similar  to  that  leading  to  Fig. 13,  it  was  found  that  a 

percentage  error  of  2X  in  r  provided  the  best  correction  at  the  first  peak  in 

tan  An.  It  is  remarkable  that  the  shape  of  the  correction  curve  at  higher 

s 

frequencies  was  very  similar  to  that  computed  from  an  assumed  error  in  R; 
thus,  experimentally  one  cannot  decide  between  a  constant  error  in  r  or  in  R 
as  the  cause  of  the  structure.  Of  course,  a  better  than  2X  accuracy  in  the 
measurement  of  the  masses  of  tip  and  sample  is  easily  obtained,  which  makes  an 
error  in  r  a  less  probable  cause  for  the  observed  structure  than  an  error  in 
R. 


A  relative  error  in  the  measured  length  of  the  bar  translates  directly 
into  an  equal  relative  error  in  a  »  k£.,  without  any  amplification. 

An  error  in  the  ratio  R  may  have  various  causes.  First,  it  may  be  due  to 
an  error  in  the  calibration  of  the  accelerometers;  even  introducing  the 
corrections  to  the  nominal  sensitivities,  as  a  function  of  selected 
frequencies  (provided  by  the  manufacturer),  may  not  be  sufficiently  precise. 
In  the  second  place,  it  is  quite  possible  that  the  lack  of  rigidity  in  the 
attachment  of  the  bar  to  tip  and  base  could  cause  an  error,  which,  when 
amplified  by  the  above  effect,  adds  to  the  observed  structure. 

To  investigate  this  source  of  error,  a  model  was  developed  for  the  wax 
connections  at  the  base  and  tip  of  the  bar,  described  in  Appendix  B.  With 
thicknesses  of  1  mm  of  wax  at  both  tip  and  bar  or  either  one  separately, 
errors  in  R  were  derived  of  about  the  right  order  of  magnitude.  In  Figs.  16 
and  17  the  value  of  the  relative  error  in  R  is  shown  as  caused  by  a  1  mm  layer 
of  wax  at  the  tip  and  base,  or  on  either  surface  separately.  The  values  for 
the  dynamic  Young’s  modulus  as  measured  at  NRL-ISRD  were:  E's  1.39  GPa, 
tan  ^£^0.015,  and  ^1520  kg/m^. 
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ncoucNcr  (Mt) 


rig.  16.  MagnitttdM  of  tka  ralatir*  arrer  ia  tba  ratio  B,  lAt/ll,  eo^otad  from  aedal  for 
wax  eonaaetioao  at  tka  kaoa  (tkiekaaoo  d^)  aad  tip  (tkiekaaoo  d^)  of  tka  saapla  kar . 
Solid  carrai  d^  a  1  am,  d^  «  1  aat  daokad  earra,  *  "**  dottad-daokad 
earraf  d^  ■  0  aat  d^  *  1  aa. 


*  -90.0 


FRtOUtHCt  <kHl) 


rig.  17.  rkaoa  of  tka  ralatira  arror  ia  tka  ratio  B|  ^ (Ab/B)  t eeapatad  froa  aodal  for  wax 
eoaaaetioao  at  tka  kaoa  (tkiekaaoo  d^)  aad  tip  (tkiekaaoo  d^)  of  tka  oaapla  kar.  Solid 
earva,  d,  ■  1  aa.  d^  >  1  aa|  daokad  earro,  d.  *  1  aa,  d^  ■  6  aa|  dottad-daokad  earra,  d 
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The  absolute  value  of  the  relative  error  due  to  the  mx  reaches  values  up 
to  15X.  It  is  remarkable,  though,  that  the  curves  for  the  relative  error  in  E' 
(Fig.  18)  and  the  error  in  tan  6  (Fig.  19),  found  by  applying  the  correction 
due  to  the  wax,  are  very  similar  to  those  for  constant  AR/R  or  Ar/r.  They  show 
that  the  corrections  for  E'  and  tan  6^  at  the  first  structure  are  of  the  right 
order  of  magnitude,  but  they  are  not  sufficient  to  explain  the  higher 
frequency  peaks. 

The  conclusion  to  be  drawn  from  this  error  analysis  is  that  the  method  of 
measuring  the  dynamic  Young’s  modulus  is  suitable  for  reaching  considerable 
precision,  on  the  order  of  one  percent,  even  in  the  frequency  range  well  away 
from  the  antiresonances.  In  the  first  place,  it  is  necessary  that  the 
sensitivity  of  the  accelerometers  is  sufficiently  well  known.  It  would  be 
possible,  in  principle,  to  eliminate  the  effect  of  inaccuracy  in  the 
sensitivity  of  the  accelerometers  by  measuring  the  ratio  R^^  with 
accelerometers  1  and  2  (sensitivities  s^  And  S2)  attached  to  tip  and  base 
respectively,  and  then  finding  R2  after  the  two  accelerometers  were  switched. 
Given  the  accelerations  a^^  and  a2  of  tip  and  base,  one  finds  that  R^^  = 

a^.  The  desired  ratio  of  accelerations  of  base  and 
tip  would  be  given  by  R  »  V(Rj^R2)>  This  has  not  been  tried  in  practice.  In  the 
csecond  place,  it  may  be  necessary  to  account  for  the  wax  connections  at  tip 
and  base.  It  is  possible,  of  course,  that  other  sources  of  error  exist  that 
mimic  the  effect  of  the  error  sources  discussed  above.  (One  might  think  of 
effects  due  to  the  discontinuities  in  the  bar  at  the  tip  and  base,  resulting 
in  effects  not  described  in  the  one-dimensional  theory.)  Moreover,  errors  from 
more  than  one  source  may  be  present  simultaneously. 
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It  would  seem  to  be  a  valid  interpolation  to  correct  the  observed  ratio  R 
by  a  constant  percentage  error  for  each  of  the  ranges  between  two  successive 
antiresonances.  This  interpolation  has  not  been  performed  in  the  present 
study,  except  for  the  above  examples.  Since  the  mun  purpose  of  this  study  was 
to  examine  the  behavior  of  various  combinations  of  the  constituents  in  a 
composite  alumer,  the  representative  values  of  E'  and  tan  6^  have  been  taken 
at  a  fixed  frequency  of  708  Ha  (this  frequency  was  a  measurement  point) ,  where 
the  deviation  due  to  uncertainty  in  R  is  small.  In  general,  this  is  close  to 
the  location  of  the  first  antiresonance.  The  frequency  of  the  first 
antiresonance  is  given  in  Tables  4a  through  4d. 

In  Ref.  10  it  was  reported  that  the  Young’s  modulus  measurement  as  a 
function  of  frequency  for  foamed  aluminum  showed  a  negative  trend.  This  slope 
could  be  reduced  to  sero  by  introducing  a  small  percentage  shift  in  R.  Since 
the  velocities  of  tip  and  base  were  measured  by  laser  Doppler,  the  error 
cannot  be  due  to  accelerometers  in  this  case. 

The  choice  of  the  ratio  r  »  ^^^^2  mainly  determined  by  the  lower 
frequency  limit  that  one  wants  to  attain.  For  instance,  with  the  values  for 
the  parameters  given  above  and  a  1%  error  in  R,  a  IX  error  in  E'  will  be  found 
at  a  frequency  of  407  Hz  with  the  ratio  r  »  1.  For  a  ratio  of  r  -  0.1,  this 
one  percent  error  in  E'  is  found  at  a  frequency  of  665  Hz;  for  r  =  10,  the 
frequency  limit  for  errors  in  E'  below  IX  decreases  to  149  Hz. 


RESULTS  FOR  YOUNG’S  MODULUS 

In  Table  4  the  results  for  the  real  part  of  Young’s  modulus  E'  and  its 
loss  tangent  tan  6^  are  collected.  Table  4  has  four  parts  (a,  b,  c,  and  d)  for 
the  four  polyurethanes  PR  1526,  PR  1538,  PR  1574,  and  PR  1590,  respectively. 
In  each  case  the  data  for  the  pure  polymer  are  given  first,  computed  from  the 
polynomial  coefficients  in  Ref.  7,  followed  by  the  Young’s  moduli  for  the 
foamed  aluminum  taken  from  Ref.  10. 
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Table  4  -  Values  of  Young’s  modulus  for  polyurethanes,  foamed  aluminum  and 
alumers.  The  ralues  for  the  real  part  E'  and  loss  tangent  tan  6^  are  giren 
for  a  frequency  of  708  Hs.  The  lowest  antiresonance  frequency  f^  is  stated 
where  applicable,  to  an  accuracy  determined  by  the  sampling  density.  Unless 
otherwise  indicated,  the  temperature  is  about  22*0.  The  data  for  the 
polyurethanes  are  from  Ref.  0;  the  foamed  aluminum  data  from  Ref.  6. 

There  are  four  parts  on  separate  pages:  4a  -  PR  1526;  4b  -  PR  1538; 

4c  -  PR  1574;  4d  -  PR  1580. 

Legend  (for  all  four  parts) : 


Tip 


nominal  cross  section  mass  of  tip, 
of  bar  and  tip  (incl.  accelerom.) 


s  - 

Short  steel 

3/4X3/4  in^ 

19.7 

g 

T  - 

Thick  steel 

3/4X3/4  ini 

37.2 

g 

A  - 

Short  alum. 

3/4X3/4  in^ 

8.7 

g 

SS  - 

Short  steel 

9/16x9/16  in^ 

11.8 

g 

ST  - 

Thick  steel 

9/16x9/16  in^ 
9/16x9/16  in^ 

21.4 

g 

SA  - 

Short  alum. 

5.5 

g 

Standard  length  bar  is  about  11.3  cm  (4.5  in.) 
Long  -  Long  bar,  about  14.9  cm  (5  7/8  in.) 


4a 

Conditions 

Date  E' (MPa)  tan  dg  f ^  (Hz) 

PR  1526,  p=1000  kg/m^. 

Foamed  aluminum,  6X,  10  ppi 

Alumer 

16.7  0.60 

207 

4/24/91  224.6  0.094  562 

4/29/91  212.2  0.076  562 

Average  219  0.085 
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4a  (Continued) 

Conditions 

Date 

E,  (UPa) 

tan  dg 

^ar 

Foaaed  aluainua,  6X,  20  ppi 

104  or  157 

Aluffler 

4/19/91 

288.4 

0.061 

631 

4/22/91 

301.8 

0.079 

699 

7/02/91 

284.8 

0.066 

643 

Average 

202 

0.060 

Foaaed  aluninun,  6X,  40  ppi 

174 

Aluffler 

4/12/91 

244.7 

0.112 

597 

4/10/91 

237.7 

0.106 

562 

5/01/91 

226.0 

0.008 

562 

Average 

238 

0.11 

Foamed  aluminum,  12%,  10  ppi 

236 

Alumer;  A 

8/06/91 

260.2 

0.082 

926 

S 

8/01/01 

273.2 

0.080 

794 

S 

8/02/91 

274.6 

0.081 

794 

T 

7/31/91 

269.0 

0.082 

656 

Average 

272 

0.081 

Foamed  aluminum,  12%,  20  ppi 

311 

Alumer;  A 

8/08/01 

329.2 

0.065 

1020 

S 

8/12/91 

331.0 

0.066 

858 

T 

8/13/91 

334.7 

0.070 

708 

Average 

332 

0.067 

Foamed  aluminum,  12%,  40  ppi 

251 

Alumer;  A 

8/28/91 

462.7 

0.073 

1166 

S 

8/14/91 

477.9 

0.075 

1000 

T 

8/13/91 

484.4 

0.073 

825 

Average 

475 

0.074 

S;  19*C 

8/21/91 

409.6 

0.085 

1020 

S;  21*C 

8/19/91 

488.7 

0.079 

1000 

S;  27*C 

8/20/91 

465.2 

0.067 

981 
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4b 


Conditions 

Date 

E'  (MPa) 

tan  dg 

^ar 

PR  1538,  />*1120  kg/m^, 

43.7 

0.29 

Foamed  aluminum,  6X,  10  ppi 

207 

A1  tuner 

1/25/91 

172.3 

0.143 

502 

2/06/91 

173.0 

0.133 

532 

2/08/91 

183.7 

0.162 

562 

2171191 

181.3 

0.123 

532 

7/19/91 

187.7 

0.120 

532 

A 

7/25/91 

189.2 

0.122 

810 

S 

7/24/91 

189.0 

0.122 

668 

T 

7/30/91 

191.4 

0.128 

552 

Long,  A 

7/24/91 

202.5 

0.113 

573 

Long,  S 

7/23/91 

206.6 

0.109 

502 

SA 

9/10/91 

176.9 

0.150 

750 

SS 

9/09/91 

174.9 

0.139 

631 

ST 

9/08/91 

178.2 

0.142 

501 

Long,  SA 

10/01/91 

173.9 

0.148 

573 

Long,  SS 

10/03/91 

165.9 

0.150 

492 

Long,  ST 

10.04/91 

165.4 

0.157 

414 

Average 

181 

0.136 

Long,  S 

7/17/91 

233.8 

0.103 

447 

Long,  S 

7/17/91 

231.7 

0.099 

438 

SA,  21*C 

9/18/91 

177.8 

0.151 

750 

SA,  26*C 

9/18/91 

163.3 

0.137 

722 

Foamed  aluminum,  6%,  20  ppi 

194  or  157 

Alumer 

1/31/91 

295.7 

0.071 

667 

2/04/91 

300.9 

0.088 

708 

2/08/91 

302.4 

0.082 

708 

2/27/91 

297.5 

0.079 

667 

Average 

299 

0.080 

SA,  23»C 

12/04/91 

241.9 

0.075 

858 

SS,  23*C 

12/06/91 

252.6 

0.084 

736 

Foamed  aluminum,  40  ppi 

174 

Alumer 

1/28/91 

234.2 

0.109 

631 

1/30/91 

230.0 

0.084 

631 

1/30/91 

238.0 

0.096 

631 

2/08/91 

225.9 

0.093 

708 

2/25/91 

227.3 

0.105 

708 

10/25/91 

236.6 

0.102 

750 

SA 

10/25/91 

221.8 

0.106 

825 

SS 

10/23/91 

222.4 

0.106 

695 

ST 

10/28/91 

219.0 

0.106 

584 

Average 

228 

0.100 

ST.  20*C 

11/25/91 

233.0 

0.108 

607 
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4c 

Conditions 

Date 

E‘  OfPa) 

tan  5^ 

(H») 

PR  1574,  />=1(X)0  kg/m^ 

203 

0.12 

Foamed  aluminum, 

6X,  10  ppi 

207 

Alumer 

5/00/01 

563.5 

0.047 

000 

5/16/01 

564.6 

0.046 

801 

6/03/01 

550.6 

0.056 

801 

6/17/01 

580.5 

0.050 

801 

6/20/01 

586.5 

0.058 

801 

Average 

571 

0.051 

SS,  23*C 

12/20/01 

472.4 

0.052 

1030 

Foamed  aluminum. 

6X,  20  ppi 

104  or  157 

Alumer 

5/06/01 

666.0 

0.032 

1000. 

5/31/01 

633.0 

0.040 

1000. 

Average 

650 

0.041 

Foamed  aluminum. 

6X,  40  ppi 

174 

Alumer 

5/05/01 

566.0 

0.047 

801 

6/03/01 

555.0 

0.055 

801 

Average 

561 

0.051 
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4d 

Conditions 

Date 

S'  (UPa) 

tan  6^ 

^ar 

PR  1590,  />=1100  kg/a^ 

42.4 

0.41 

Foaaed  nluainun,  6X,  10  ppi 

207 

Aluaer 

3/01/91 

215.8 

0.106 

562 

3/20/91 

230.4 

0.110 

562 

3/29/91 

233.4 

0.119 

562 

4/01/91 

226.3 

0.101 

562 

6/24/91 

227.6 

0.103 

584 

Arerage 

227 

0.108 

SA  21.6*C 

11/07/91 

203.7 

0.132 

779 

ST  21.5*C 

11/07/91 

203.3 

0.132 

562 

Arerage 

199 

0.140 

SS  22.0*C 

11/01/91 

190.3 

0.167 

631 

Ponaed  aluainua,  6X,  20  ppi 

194  or  167 

Aluaer 

3/01/91 

316.4 

0.082 

708 

3/18/91 

318.7 

0.079 

708 

SA  21.5*C 

11/05/91 

322.4 

0.099 

962 

SS  21.5«C 

11/04/91 

319.0 

0.096 

810 

ST  21.5*C 

11/05/91 

319.8 

0.108 

681 

Arerage 

319 

0.093 

Foaaed  aluainua,  6X,  40  ppi 

174 

Aluaer 

3/01/91 

258.0 

0.113 

631 

3/18/91 

248.2 

0.120 

631 

ST  22*C 

10/29/91 

245.4 

0.126 

610 

Average 

250 

0.120 

SA  21.S*C 

11/08/91 

275.3 

0.122 

910 

SS  21.5*C 

11/12/91 

283.2 

0.128 

779 

ST  21.5»C 

11/13/91 

279.4 

0.127 

644 

Average 

279 

0.126 
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The  legend  for  Table  4  shows  the  warious  values  of  the  paraneters  used  in 
the  experiments.  Most  of  the  measurements  are  done  for  a  standard  bar,  with 
nominal  length  of  11.4  cm  (4.5  in.);  whenever  a  long  bar  was  used  this  is 
indicated  by  *long,'  its  nominal  length  is  14.0  cm  (5  7/8  in.).  The  parameters 
were  varied  in  order  to  check  the  validity  of  the  experiments:  the  values  for 
E'  and  tan  6^  should  be  independent  of  the  choice  of  parameters.  This  is  true 
in  most  situations  and  is  discussed  further  below. 

By  scanning  the  values  in  Table  4,  one  may  obtain  a  feeling  for  the 
spread  in  the  data.  The  values  for  Young *s  modulus  are  given  for  the  same 
frequency  (708  Hs) ,  where  the  structure  referred  to  in  an  earlier  section  is 
of  little  influence.  For  general  information,  the  lowest  frequency  is  given 
where  the  absolute  value  of  the  tip-to-base  voltage  ratio  is  maximum.  This  is 
close  to  the  first  antiresonance  of  the  bar-tip  system.  The  accuracy  of  this 
value  is  dependent  on  the  sampling  density  used. 

As  the  simplest  of  models  to  combine  moduli  of  constituents,  one  may 
consider  the  two  extremes  of  a  series  or  a  parallel  combination.  Assuming  that 
in  the  altimer  the  displacement  of  the  two  constituent  solids  is  the  same, 
while  the  applied  force  is  distributed,  one  could  opt  for  an  addition  of  the 
Young’s  moduli  of  the  constituents  to  obtain  the  modulus  of  the  composite,  for 
a  first  crude  insight  into  the  measurement  results. 

The  following  features  in  the  data  deserve  attention: 

1.  Table  4a.  For  the  alumers  based  on  polyurethane  PR  1526,  the  6X,  10  ppi 
type  has  an  E'  that  is  about  equal  to  the  sum  of  the  moduli  for  the  foamed 
aluminum  and  the  pure  polymer.  The  value  is  slightly  larger  for  the  40  ppi 
type,  but  the  value  of  E'  is  considerably  larger  than  this  sum  for  the  20  ppi 
type.  This  points  out  that  the  pore  sise  of  the  aluminum  matrix  has  a 
noticeable  influence  in  addition  to  its  porosity,  contrary  to  the  assumption 
in  the  various  models,  which  ignore  the  pore  size  effect.  The  influence  of 
pore  size  is  confirmed  in  all  the  other  data  presented.  This  point  is  brought 
to  the  fore  by  the  12X  relative  density  foam:  here  the  20  ppi  type  shows 
little  difference  in  the  E'  value  between  foam  and  alumer,  while  the  40  ppi 
type  displays  considerable  difference.  Notice  that  the  pore  size  of  the  foamed 
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aluminum  depends  on  the  relative  density,  in  addition  to  the  pores-per-inch 
designation;  for  a  higher  density  foam,  the  mall  thickness  is  larger  and  thus 
the  pore  sise  is  smaller.  The  loss  tangent  does  not  display  any  sisable 
differences  from  one  type  of  alumer  to  the  next,  but  it  is  much  smaller  than 
the  loss  tangent  of  the  pure  elastomer,  which  is  near  0.6. 

The  spread  in  the  results  for  different  tip  masses  is  well  within  the 
error  limits.  The  measurements  at  different  temperatures  for  the  12%,  40  ppi 
type  show  a  decrease  for  increasing  temperature,  as  would  be  expected  for  a 
viscoelastic  material;  the  steps  appear  larger,  though,  than  those  for  the 
pure  polymer  at  the  same  temperatures. 

2.  Table  4b.  With  alumers  based  on  the  PR  1538  elastomer,  this  table  shows  a 
reasonable  agreement  between  the  results  for  the  6W,  10  ppi  type.  Only  the  two 
measurements  for  a  long  bar  with  regular  cross  section  (Long,  S)  are  somewhat 
outside  the  expected  range  of  error,  for  unknown  reasons.  The  20  ppi  type 
again  has  an  E'  that  is  considerably  larger  than  for  the  other  two  types, 

10  ppi  and  40  ppi.  The  temperature  variation  is  as  expected.  The  consistency 
of  the  40  ppi  data  lor  various  tip  masses  at  small  cross  section  (SA,  SS,  ST) 
is  good.  The  loss  tangent  for  the  composite  is  again  considerably  smaller  th2m 
that  of  the  pure  elastomer. 

3.  Table  4c.  The  polymer  PR  1574  is  much  sore  rigid  than  the  other  three;  this 
apparently  leads  to  a  large  modulus  for  the  alumer,  larger  even  than  the  sum 
of  the  moduli  of  the  constituents.  The  20  ppi  type  again  stands  out  by  its 
larger  E'  as  compared  with  the  10  ppi  and  40  ppi  types.  The  loss  tangent  shows 
little  variation  between  the  alumers;  its  value  is  less  than  one  half  of  that 
of  the  pure  polyurethane. 

4.  Table  4d.  For  the  alumers  based  on  PR  1590,  one  finds  again  that  the  20  ppi 
type  has  a  larger  E'  than  the  other  types.  The  loss  tangent  is  three  to  four 
times  smaller  than  that  of  the  pure  elastomer.  For  the  SA,  SS  and  ST 
measurements  of  the  10  ppi  type,  at  various  temperatures,  the  values  are 
contrax7  to  what  follows  from  viscoelastic  behavior:  the  value  at  21.5*0  is 
smaller  than  at  22*0,  although  not  by  much.  The  cause  of  these  deviations  is 
unknown.  The  variations  due  to  temperature  in  the  20  ppi  type  are  w.'thin  the 
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limits  of  error;  for  the  40  ppi  type  the  effect  of  0.5*  difference  in 
temperature  is  quite  noticeable  and  has  the  correct  sign;  the  lower 
temperature  results  in  a  larger  E' . 

DlLmnOHAL  1AVB  SPEED  IN  ALONERS 

Table  5  shows  rarious  properties  of  alumers  and  the  polyurethanes  on  which 
they  are  based.  The  bulk  modulus  is  taken  from  Table  1,  with  interpolation  of 
the  temperature  wariation  to  approximate  the  ralue  at  22*C.  Young’s  modulus  for 
the  pure  polyurethanes  is  computed  from  the  polynomial  coefficients  in  Ref.  7. 
Young’s  modulus  for  the  alumers  is  taken  from  Table  4. 

For  application  as  an  acoustically  transparent  material,  the  dilatational 
wave  speed  c^,  in  addition  to  the  density  /),  is  important.  The  lowest  ware 
speed  is  found  for  the  PR  1538,  6X,  20  ppi  alumer  (Table  5).  Its  density  is 
greater  than  that  of  the  PR  1526  types,  and  thus  the  three  types  of  alumer 
based  on  the  PR  1526  polymer  hare  the  lowest  ^c-product,  close  to  1.7  SI  Mrayl. 
This  ralue  is  only  3X  smaller  than  the  ^c-walues  for  the  PR  1538  alumers.  The 
non-dimensional  loss  coefficient  a/k'  raries  by  about  a  factor  of  three  in  the 
set  of  aliiners;  it  is  smallest,  0.0033,  for  the  PR  1574,  6X,  20  ppi  type  and 
largest,  0.010,  for  the  PR  1560,  6X,  40  ppi  type. 
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Table  5.  Various  computed  elastic  properties  of  polyurethanes  and  aluaers  at 
20*C. 

/^-density;  6' -real  part  of  shear  modulus;  i^',i/"-real  and  imaginary 
part  of  Poisson’s  ratio;  c^,  s/k'-vaTe  speed  and  non-dimensional 
loss  coefficient  of  dilatational  mayes;  the  complex  ware  number  is 
k*  »  k'  -  is;  ~  product  relatire  to  seawater  with 

pc  ^  1.539  SI  Urayl. 


Material 

P 

kg/m3 

G' 

MPa 

1/' 

m/s 

PCj 

rel 

Polymer 

PR 

1526 

1000 

5.57 

0.4988 

0.0014 

1525 

0.0059 

1.525 

0.991 

Alumer 

10 

ppi 

1068 

73.7 

0.4814 

0.0020 

1584 

0.0088 

1.691 

1.099 

6X 

20 

ppi 

1070 

98.5 

0.4814 

0.0022 

1608 

0.0078 

1.718 

1.116 

40 

ppi 

1057 

79.4 

0.4851 

0.0031 

1702 

0.0069 

1.702 

1.106 

Polymer 

PR 

1538 

1120 

14.6 

0.4969 

0.0018 

1455. 

0.0042 

1.630 

1.059 

Alumer 

10 

ppi 

1121 

60.8 

0.4886 

0.0029 

1558 

0.0069 

1.746 

1.135 

6X 

20 

ppi 

1152 

101, 

0.4804 

0.0028 

1526 

0.0097 

1.758 

1.142 

40 

ppi 

1134 

76.7 

0.4853 

0.0029 

1540 

0.0078 

1.746 

1.135 

Polymer 

PR 

1574 

1000 

68.3 

0.4866 

0.0031 

1616 

0.0055 

1.616 

1.050 

Alumer 

10 

ppi 

1093 

195. 

0.4673 

0.0033 

1704 

0.0040 

1.863 

1.211 

6X 

20 

ppi 

1098 

222. 

0.4631 

0.0030 

1716 

0.0033 

1.884 

1.224 

40 

ppi 

1086 

191. 

0.4683 

0.0028 

1718 

0.0062 

1.866 

1.212 

Polymer 

PR 

1590 

1100 

14.2 

0.4971 

0.0023 

1497 

0.0065 

1.647 

1.070 

Alumer 

10 

ppi 

1169 

76.4 

0.4859 

0.0027 

1543 

0.0092 

1.804 

1.172 

6X 

20 

ppi 

1160 

108. 

0.4805 

0.0032 

1572 

0.0085 

1.824 

1.185 

40 

ppi 

1180 

81.2 

0.4846 

0.0033 

1548 

0.0101 

1.826 

1.186 
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CONCLUSIONS 

•The  bulk  modulus  of  m  giTen  alumer  is  larger  than  the  bulk  modulus  of  the 
polyurethane  used  as  a  filler  bj  slightly  more  than  the  relatire  density 
of  the  foamed  aluminum  matrix.  (This  rule  is  obserred  for  the  6X  aluminum 
type  only.) 

•The  data  set  is  insufficiently  large  and  the  available  models  inadequate  to 
reliably  predict  the  Young’s  modulus  of  new  combinations.  Systematic 
effects  were  found  that  point  to  the  importance  of  the  pore  sise  in 
determining  Young’s  modulus  of  the  alumer. 

•The  loss  tangent  of  Young’s  modulus  of  an  alumer  is  generally  much  smaller 
than  the  loss  tangent  of  the  constituent  pure  polymer. 

•The  pc-values  for  the  alumers  studied  vary  from  1.60  to  1.87  SI  Mrayl. 

•A  preliminary  study  of  possible  sources  of  error  in  the  measurement  of  the 
dynamic  Young’s  modulus  suggests  that  a  high  precision  may  be  attained. 
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Appendix  A 

STIFFNESS  AND  COMPLIANCE  MATRICES 


The  following  conpliance  and  stiffness  satrix  relations  applj  to  a  system 
with  one  rotational  symmetry  axis  [15] ,  which  is  identical  to  the  hexagonal 
system.  The  symmetry  axis  is  indicated  by  the  subscript  3. 

The  elastic  conpliance  matrix  is: 


“11 

“12 

“13 

0 

0 

0 

“12 

“ll 

“13 

0 

0 

0 

“l3 

“13 

“33 

0 
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“44 

0 

0 

0 

0 

0 

0 

“44 

0 

0 

0 

0 

0 

0 

“44. 

(Al) 


®66  ”  ^ 

The  stiffness  matrix  [c]  has  the  same  fora  as  the  matrix  [s] ,  but  here 
Cgg  =  0.5 

Relations  between  elements  of  the  stiffness  matrix  [s]  and  the  compliance 
matrix  [s]  are 

“33 

°11  *  ®12  *  T 


®11  ■  <=12  = 


*11"“12 


c,-  = 


'13 

<=33 
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“ll^“l2 


c..  =  s 


44  "44' 


with  s  «  Sgg  (■ii**i2)  "  2  *13^ • 


(A2) 
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Appendix  B 

MATHEMATICAL  MODEL  FOR  EFFECT  OF 
WAX  CONNECTIONS 


In  this  Appendix  corrections  are  derived  for  the  ratio  of  accelerations  of 
base  and  tip  of  the  sample  bar,  due  to  the  lack  of  rigidity  of  the  max 
connections . 


sample 

< - <*2 - ^ 

max  max 


Figure  B1  shoms  the  basic  featiires  and  definitions  of  some  parameters  of 
the  model.  A  bar  sample  is  attached  by  tmo  max  layers  to  the  tip  and  base. 
From  linear  maveguide  theory  one  obtains  the  equations  for  the  forces  per  unit 
area  f^j  and  velocities  v^.  at  the  interfaces  in  the  fora 


^11  ^11  '^11  ^12  ^21 
^21  =  ^12  ^11  *  ^11  ^21 
^12  “  *11  ^12  *  *12  ^22 
^22  ^  *12  ^12  *  *11  ^22 
*13  “  *11  *13  *  *12  *23 
*23  “  *12  *13  *  *11  ^23* 
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^11  *  i  ) 

^12  =  ^I'^l  /  (  i  > 

*11  =  ^2*^2  /  (  i  tan  kjdg  ) 

*12  *  /’2®2  /  (  ^  ^ 

*11  “  ^1®1  /  (  i  tan  kjdjj  ) 

*12  =  ^1*"!  /  C  i  ) 

d^  -  thickness  of  nax  at  tip 
■’  length  of  sample  bar 
In  Eq.  (B2)  one  has,  for  a  harmonic  ware. 


Pi 

P2 

*=2 

>^1 

dh 

“2 


density  of  wax 
bar  speed  in  wax 
density  of  sample 
bar  speed  in  sample 
warenumber  in  wax 
mrenumber  in  bar 
thickness  of  wax  at  base 
mass  of  the  sample  bar 


-  ^  “3  ^21  .  ^ 

^21  "  A  ~  ‘*3  ^21’ 


(B7) 


where  w  is  the  angular  frequency,  m^  the  mass  of  the  tip,  and  A  the  cross* 
sectional  area.  The  forces  and  relocities  at  the  two  interfaces  of  the 
interfaces  of  the  sample  bar  are  simply  related,  and  thus  the  original 
unknowns,  six  forces  and  six  relocities,  reduce  to  eight  unknowns.  Eq. (B7) 
reduces  the  number  to  sewen.  Equation  (B5)  is  extraneous  to  the  computation, 
and  thus  five  homogeneous  equations  remain  in  six  uxiknowns. 


One  may  solve  Eqs.  (Bl,  B2,  B3,  B4,  and  B6)  for,  say,  fj^^,  fj^2>  ^13* 

and  V2j^  in  terms  of  y^2' 

The  observed  ratio  =  v  (base) /v  (tip)  =  -  ~  ^”^12/^21^  ^^13/^12^* 

R  is  related  to  the  ratio  used  in  the  data  reduction  R  byR  =  R  +  AR, 
o  c  '  c  o 

where 


*11‘*‘*3 

“c*  k2'*2  - 


From  the  algebraic  solution  to  the  equations  it  follows  that 

-V, 


112 

^21 


^  2  ^  ^c^yii~yi2^  ^  ^11  "^12  ^ 
®  712  712  *12 


*11  *3  1 

“12J  U12J 


R-+ 


(B8) 


(B9) 


Retaining  terms  of  order  and  k^d^,  qj^  is  approximated  by 
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2  _  2 


.111212  f!u]  f!3.]. 

^  ^12*12  KaJ  {^12} 


(BIO) 


Next,  the  solution  gires  the  second  factor  in  the  fon 


^13  _  ^11  ^ 

’12  *12  * 


(Bll) 


uhere  the  first  tern  on  the  right-hand  side  is  of  order  1,  and  A  is  shorthand 
for  a  correction  tern,  to  be  stated  explicitly  later. 

Thus  one  aay  urite: 


».  =  (S  *  *1  ■  “i*  fe  ■ ') 

+  A  q^. 

The  correction  ten  q2  is  found  froa  the  algebn  to  be: 

*12^12  *12^12 


Retaining  only  errors  of  order  k^d^  and  k^d^^,  q2  is  approxiaated  by 


2  _  2 


A  *  f  \  f  * 

2  11~*12  *11  *3 

^  *12  *12  [*i2J  riaJ  * 


(B12) 


(B13) 


(B14) 


The  correction  ten  q^^  is  siailar  in  fon  to  the  correction  ten  q2i  they 
represent  errors  due  to  the  wax  at  the  tip  and  base,  respectirely. 
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SDIFLIPIED  EXPBESSIONS  FOE  GORBBGTIONS 


In  the  following  deriyations  it  will  be  assumed  that  uid  h^d^  are  auch 

smaller  than  one,  so  that  sin  hj^d^  s  ^l^t  ~ 

By  inserting  the  definitions  for  the  quantities  y..,  and  z,,  the  first 

xj  13  o 

term  in  the  expression  for  may  be  reduced  to  a  simpler  form  as  follows: 


2  2 
^12  *12 


Pi 

1^2J 


Vt 


(BIS) 


and  the  second  term  leads  to 


ul  ffa.! ,  (-3/*) » fci-it 

I2J  [j'lzJ  I 


cos 


!3 

L“U 


(kid^)-*  cos  k2d2 


(B16) 


f"3l  P2l 

=  ra  TO  ^2*^2 

where  E2  are  the  Young’s  moduli  of  the  wax  and  the  sample,  respectiyely . 


In  a  similar  fashion  one  may  reduce  the  first  term  of  the  expression  for  ^ 
as  follows: 


.2  2 

*ir*i2 

*12*12 


^2°2 


l^l^lJ 


fEa) 


kidbSin  kgdj  = 


P2^2 


2]  ,  [f2 

I)  ^  1*^1. 


1®1J 


(B17) 


Vb**“  V2’ 

and  the  second  term  of  q2  becomes 


mu.  MmoiAHSoii  upoit  tsm 


The  total  correction  AR  is  equal  to  thus 

“  =  ^  ^  “•  ^dj,«  ^  dj  .In  kjdj. 


(Bl«) 
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